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Abstract 

It is demonstrated that due to back-reaction of quantum effects, ex- 
pansion of the universe stops at its maximum and takes a turn around. 
Later on, it contracts to a very small size in finite future time. This 
phenomenon is followed by a "bounce" with re-birth of an exponentially 
expanding non-singular universe. 

PACS no. 98.80.Cq 
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Astrophysical data [1, 2] show cosmic acceleration in the current uni- 
verse, which is fuelled by dark energy (DE) violating strong or weak 
energy condition (SEC or WEC). If DE mimics quintessence, the equa- 
tion of state parameter (EOSP) w > — 1 and SEC is violated . In the 
case of phantom (super-quintessence) DE, w < — 1 and WEC is vio- 
lated. It is found that when w< —1, DE leads the cosmological model 
having "big- rip singularity" and ending in cosmic doomsday with infi- 
nite energy density [3] . A comprehensive review on DE and accelerating 
universe is available in [4] . In contrast to quintessence model, phantom 
model shows growing curvature with time and in models with "big-rip 
singularity" at finite time ts with t being the cosmic time, curvature in- 
variants (containing high powers of curvature) grow very strong, when 
|is — i| is sufficiently small. Quantum gravity suggests that quantum 
corrections to matter fields like scalars, spinors and vectors, in curved 
space-time, depend upon curvature. It implies dominance of quan- 
tum gavity near t = ts- It is analogous to the early universe model, 
where curvature invariants are very strong and cause effective role of 
quantum corrections. This fact was pointed out in [5, 6], where it is 
demonstrated that singularity may be avoided or made milder using 
quantum corrections near future singularity time ts- 

Like [6], here also, scale factor of the universe is assumed in the 
beginning. Further, it is probed what type of conformal scalar (quin- 
tessence or phantom) support this scale factor (the assumed one) for 
being the solution of Priedmann equations. Like [5, 6], here also, quan- 
tum gravity effects are investigated when t approaches ts- In [6], a 
possible escape from cosmic doomsday at t = tg is discussed using 
quantum corrections. Here, the approach is more elegant and a cosmic 
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scenario is obtained, where expansion stops at t = tm very near to tg. 
As a result, it takes a turn around and contracts for a small period 
{ts — tm)- Subsequently, at t — ts, cosmic bounce takes place and uni- 
verse expands exponentially for t > tg- This scenario emerges due to 
quantum effects during the period {ts — tm)- Thus, here, singularity is 
not made milder , but it is completely avoided in contrast to results 
in [5, 6]. In the present paper, this scenario is obtained for (i) general 
relativity (GR) model, (ii) Randall- Sundrum II (RSII) model and (iii) 
Gauuss-Bonnet(GB) model. 

It is found that the scale factor reaches its maximum in finite time tm, 
when scalar is phantom (super-quintessence) in GR and RSll models. 
But, in GB-model, similar situation is obtained for both quintessence 
conformal scalar and conformal phantom subject to different condi- 
tions. In [7] also, conformal phantom scalar, in FRW universe, had 
been considered in a different context. Here, an escape from 'future 
singularity' in phantom model is explored using quantum gravity ef- 
fects of conformal phantom scalars, whereas, in [7] , non-integrability of 
hamiltonian system of these fields is discussed using Ziglin theorems on 
integrability. Natural units {H — c — 1) are used with h and c having 
their usual meaning. GeV is the fundamental unit and cosmic time t 
is measured in GeV~^. 

General Relativity model of the future universe 
(a) Classical approach 

According to experimental probes [8], our universe is spatially fiat 
with growing scale factor a{t), given by the distance function 
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dS^ ^ dt^ - a\t)[dx^ + dy^ + dz^] (1) 

such that a = ^ > 0. 

DE density p and pressure p are given by Priedmann equations 

^ = U = — " 

and 

d AttG , , , , , 

- = ^(p + 3p), (26) 

where G — Mp^ [Mp is the Planck's mass). 

Here, source of the dark energy is conformal scalar (pi^x) given by the 
action 

S^^ j d'x^]^ [uj^g^^d^(t>d,(t> - ^R(t>^] , (3a) 

where u"^ = ±1 for quintessence and phantom (super-quintessence) 
scalars respectively 

This action yields the Klein-Gordon equation 

uj^4> + 3uj^H^ +{H + 2H^)(t) = 0, (36) 

where R — 6{H + 2H'^) and (f) — (t){t) due to homogeneity of the space- 
time (1). 

The energy density p and pressure p are obtined from the action (3a) 

as 



+ (3c) 



and 



P = - ^(# + 0^ + 3W) - ^(2i? + 3i/2)0l (3d) 
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The conservation equation is 

p + 3H(p + p) ^0. (3e) 

It is noted that the Klein-Gordon equation (3b) can be derived con- 
necting eqs.(3c), (3d) and (3e). So, eqs.(3b) and (3e) are not indepen- 
dent implying that it is enough to use either of these two equations. 
Now, we are left with three equations (2a), (2b) and (3b) as well as 
(3c) and (3d) for p and p respectively. Connecting (2a) and (3c), it is 
obtained that 



This equation shows that when u;^ — —1,0^ > 3/8ttG. Using 
eqs.(3b), (3c), (3d) and (4a) in eq.(2b), it is obtained that 



+ 2i?02 + X = (46) 

with 



2^=-l,[i7+(3-a;^)i/2](6u;2-5)i/2 + (u;2-l)ij]-4^u;2(L2-l)i/2 



and 



^ (3-^2)1^2^'' 



A-kG 



(4e) 



Eq.(4b) is obtained writing a/a — H + H'^as well as using 0, in terms 
of 0, if and H. This equation yields 



2 ± _ 



(4/) 
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Thus it is obtained that an a(t) will satisfy eqs.(2a) and (2b), if 
conformal scalar obeys eqs.(4a) and (4b) [9]. Now, we can use condi- 
tions (4a) and (4b) to examine an a{t) for being solutions of Friedmann 
equations (2a) and (2b) with conformal scalar as a source of DE. 

For the time period to < t < tm, the form of a{t) is assumed to be 

a{t) + B{Mpty + \Mp{ts - t)\-\ (5a) 

where r > 1 to get a{t) > 0, > 0, ^4 and B are arbitrary constants. 
The classical approach shows that the scale factor (5a) has big-rip 
singularity aX t = ts [10]. In the following section, it is shown that 
an escape from this catastrophic situation is possible using quantum 
gravity effects. 

Now, we explore whether mimics quintessence or phantom, when 

a(t) (given by eq.(5a)) satisfies cqs.(2a) and (2b) as well as acquires its 
maximum at time t = tm- It turns out to explore 0, when a(t) satisfies 
eqs.(4a) and (4b) for to < t < tm with a{tm) being the maximum. With 
Oo = a{to){to being the present time) A is obtained as 

A = ao- BiMptoY - \Mp{ts - to)\~\ (56) 
a{t), given by eq.(5a), expands to its maximum by the time t — tm, 

if 

a{tm) = vBMpt^m^^ + qMp''\ts - t^l^"^-') = 

implying 

qM'f''-'-\ts-tm\^-''-^^ 

' I'm 

Thus ,a(t) from eq.(5a), looks like 
a{t) ^ ao-^^^\ts-tmt''-''\f -fo)+M-''\ts-t\-'^-M-^\ts-to\-'' 

^ tm 

(7) 
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for to <t < tm- 
Moreover, 

a{t) = qMp\ts - tml^-"-^^ [ - {t/tmY + {\t 



(8) 

and 

a{t) = qMp% - tj-^-'^ [ -t^(^t/tmY~'^ 

(q + 1)/' \ (-9-1)1 

+ j^—^^(\ts-t\/%-t^\) J. (9) 

These equations show a{t) > and a{t) > for to < ^ < and 
5 > 0, r > 1 as 

' >i. (10) 



ts-t\ t 



Eqs.(8) and (9) yield that, at t = i7 = and 7^ 0. 

Now, wc have two cases (i) cj^ = 1 and (ii) uo"^ = —1. Eq.(4f) yields 
0^ being indeterminate for lu'^ — 1, when H — and H ^ at t — t^ . 
But, for o;^ = —1 

0^(U=3|j (11) 

for H{tm) = and H{tm) 7^ 0. 

It means that the scale factor a{t), given by eqs.(5a), (5b) and (6) 
or eq.(7), is not a solution of eqs.(2a) and (2b) for the time period 
to < t < tm, ii (p is a quintessence scalar, but it is a solution of eqs.(2a) 
and (2b) for the time period to < t < tm for (p being a phantom scalar 
characterized by a;^ = —1. 

This analysis suggests that if universe is dominated by quintessence 
dark energy, expansion will not stop at its maximum. So, in this model 
(GR model), only the case o;^ = — 1 is investigated. 
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Quantum effects and bounce near tg 



In the scenario discussed above, a{tm) — 0, which means that either 
p in eq.(2a) vanishes at t — tm or there should be a term pcorr on the 
r.h.s. of the equation (2a), which remains ineffective during the time 

interval to < t < tm, but a.tt = tm it makes the effective energy density 

PeS = P + Pcorr = 0. 

Using eqs.(4a), (5a,b) and (6) in eq.(3c), it is found that p increases 
with time during the interval to < t < tm- So, only second possibility 
may be valid as p(im) 7^ 0. This argument modifies the Priedmann 

equation (2a) as 



0'- 



(P + Pcorr) (12a) 



and eq.(2b) as 

a AttG 



{P + Pcorr + 3p + 3pcorr ) , (126) 



a 3 

such that for t < t^., Peon ~ and at t = t^, P = — Peon- In what 
follows, we explore Pcon satisfying these requirements. 

As it is found in eq.(lO) that tm also satisfies the inequality \ts—tm\ < 
1. So, when tg ^ t > tm, curvature becomes very strong and quantum 
gravity becomes effective [5, 6, 11]. Conformal anamoly for scalars, due 
to quantum effects, yield correction to p and p as 



and 



TV • TV (JH 

^ 360(47r)2^ ^ 180(47r)2^ dt ^' 

(14) 

where TV is the number of scalars. 
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Now, connecting eqs.(5a), (5b), (13) and (14), for to < ^ < ^m, it is 
obtained that 



{q + l) 



(r-1) 



\t -t I 

t 



\ {'•-1) 

t/tjj \ts-t\{q + 2) 



(15) 



This equation shows that for i < i^, 0. But t — tmi 



PA{tm) = 



NqMpllf |(-2.-4) 



{q + l) 



-{r - 1) 



\t — t 



(g + 2) 



(16) 



Thus, we find that pa satisfies the criteria required by Pcorr- So, it is 
set that Pcorr — Pa- 

As H{tm) = 0, eq.(13) yield p(tm) + pA{tm) = 0, so 

NqM~'^'^ 



p{tm) 



60(47r)2 



1/ _/ I (-25-4) 



[q + l) 



-{r - 1) 



iq + 2) 



n 2 



(17) 



Thus it is obtained that DE density grows with time . Moreover, 
the correction term which is neghgible for t < tm becomes very strong 

at t = tm- As Pa < 0, Peff(^m) = P{tm) + pA{tm) = , but for t < tm, 

Pes{tm) — P- So, it is quantum effect, which stops expansion of the 
universe at t — tm- 

For the modified Friedmann equations (12a) and (12b), solution is 
taken as 

a(t) = exp[B'{Mpt)'' + C'M^{t - tm)^ + M],\ts - tp], (18a) 

when t > tm- Moreover, a{t) should be continuous at t — tm, which is 
given by the condition \im^^od{tm + e) — \ime-^o a{tm — e)- Using this 
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condition for a(t), given by eq.(5a) (when t < tm) and a{t), given by 
eq.(18a) (when t > tm), it is obtained that 



H = aB'M^t^''-^^ + pC'M^it - tm)^^'^^ - lMl\ts - t\^^-^\ 
H = a(a-l)5'M^i("-2)+/3(/3-l)C"M^(i-i^)('^-2) 

+7(7-l)M^|i,-i|M, 
H = a{a-l){a-2)B'M^t^"-^^ + PiP-l){(3-2)C'M^{t-tm)^'^'^'^ 

-7(7-l)(7-2)M^I^.-^l^^-'\ 

H = a(a-l)(a-2)(a-3)fi'M^t("-^)+/3(/3-l)(/?-2) 

x(/3 - 3)C'M^p{t - tm)^^-'^ + 7(7 - 1)(7 - 2)(7 - 3)M^|i, - 



If a{t), given by eq.(18a), satisfies eqs.(12a) and (12b), p and p are 
obtained as 



yielding 



(186) 



For t > tm, eq.(18a) yields 



(19a, b, c, d) 




3H^ N 




SttG 180(47r) 




-[2H + 12HH + 18i/2ii' + m% 



(20a, 6) 



where H, H, H and if are given by eqs.(19a,b,c,d). 
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Eqs.(19 a,b,c,d) and (20a,b) show that p and p can be finite at t = ts, 
if 7 > 4. So, here onwards, 7 = 4 is set in equations. Moreover , 
without any harm to physics, a and (5 are also taken as a — (5 — "i. 

Thus a{t), given by eq.(18a),is obtained as 

a{t) = exp[B'{Mptf + C'Ml{t - tmf + M%\ts - t\% (21a) 

where t > t„i. If a(t), given by eq.(22a), acquires its minimum a^ at 
t^ts, a{ts) = yielding 

Moreover, eq.(18b) is re-written as 



4:Mp\ts-tr, 



|2 



B' = ^^'^"'3;, ■ (21c) 



m 



Use oi a — P — 3 and 7 = 4, in eq.(19a), yields H < implying 
a{t) < o during the period tm < t < tg- It shows contraction of the 
universe during this period. For t > tg, eq.(19a) yields 



H ^ 4M'p\ts-tmn-' 



+ 



^m' ^\ps ^m\' ^^m' ^ 



> 
(22) 



showing expansion . Ki t — tgi eq.(21a) yields 



a{ts) = exp "^M^^(t, - tmf] . (24) 



Thus, it is obtained that, during the period universe 
will contract and will acquire its minimum at t = tg. Later on, it will 
expand for t > ts showing a bounce t — ts- 

As H{ts) — a{ts)/a{ts) — 0, eqs.(13 ) and (19 a,b,c,d) yield 

P{ts)-^,MUts-tm)\^y (24a) 



12 

and 
p{ts) 



S.K.SRIVASTAVA 



JV 2M« 

607r2 ^ 



(246) 



RS II brane-world model 



Classical approach 

It is mentioned above that curvature invariants in future universe, 
dominated by phantom dark energy, get stronger with increasing time. 
So , relevance of high energy physics increases [12]. With this view, 
above type of analysis is done in brane-world also. Here, RSII model is 
considered due to its simplicity and geometric appeal. Moreover, it is 
more relevant to cosmology compared to RSI model. Modified version 
of eqs.(2a,b), for brane-gravity, are written as [6, 13] 



-PRS 



1 + 



Prs 



2A J 



(25a) 



and 



a 

a 



Prs + ^Prs + jPrs{'^Prs + ^Prs) 



(256) 



where A is the brane-tension. 
Eq.(25a) yields 



Prs 



A 



1±\ 1 + 



167rA 



(26) 



Eqs.(3c) and (26) yield 



-H4>±\ {1- uj^)H'^(l)^ + 2uj^X 



-1±A 1 + 



47rA 



(27) 
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In this model also, matter field is given by the same scalar field 0, 
so for prs and prs also eqs.(3c) and (3d) are used in eq.(25b). As a 
result, it is obtained that 



{(t-i)[-"/ 



IGttA 



X 



-1±\ 1 + 



3M|,i/2 



IGttA J 



-UJ 



3M2//2 



]] 



3+if^ + {y(if+«^)-(3+^)}[l-l±Vl + ^]]- 

(28) 

a(t), given by eq.(5a), acquires its maximum aXt — tm, so a{tm) — 0, 
but d{tm) 7^ 0. It yields H{tm) 7^ 0. In case, uo"^ = 1, eq.(28) implies 
that either H vanishes or it is complex. So, like GR-based model , a{t) 
does not satisfy modified Friedmann equations (25a,b) for RSII model 
too, when lo'^ — 1 . But, in case of a;2 = —1, a{t), given by eq.(5a) 
satisfies eqs.(25a,b) and a,t t — tm, H{tm) — but H{tm) 7^ and 

3 



(29) 



AnGil + A) 

Thus, like GR-modcl in RSII-model too, the scale factor a{t), given 
by eq.(5a) satisfy eqs.(25a,b) for a;2 = — 1 (phantom conformal scalar), 
but not for a;2 = 1 (quintessence conformal scalar). 

Quantum bounce near t — ts 

As H{tm) = 0, eq.(25a) yields 



PRs{tm) = 0, 



(30a) 
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or 

PRsitm) = -2 A. (306) 

These results are obtained without making quantum correction in eq.(25a). 
Eq. (30a) corresponds to GR- model. Effect of the brane theory is real- 
ized in eq.(30b). So, for RSII model, only eq.(30b) is considered. 

Like GR-model, here also quantum correction to eq.(25a) is made 
replacing p(^RS) by P(rs) + Pa as 

rj2 SttG r (prs + pa)] 

H = -^{PRS + pa) [1 + ^ J (31) 

Now using H{tm) — 0, in eq.(31), it is obtained that 

PRs{tm) = -2 A - pA{tm) = -2A + p{tm) ■ (32) 

where p(tm) is given by eq.(17). 

It shows that, in RSII model, expansion of the universe stops at 
a lower energy density compared to GR case due to brane tension . 
Later on,like GR case, it contracts during the interval tm < t < ts to 
its minimum a{ts) following a{t), given by eq.(21a), with B',C' given 
by eqs.(21b) and (21c). 

Like GR case, here also, a.t t = ts quantum correction to eq.(25a) 
leads to 



Ps{RS) = -2A + p{ts) (33) 

with pits), given by eq.(24a) . 

Thus, in this case too, bounce will take place att — ts like GR-model. 
In [14], bounce has been discussed using self-gravitational corrections 
to Priedmann equation for the brane cosmology. 



Gauss-Bonnet model 
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Classical approach 

Einstein-Gauss-Bonnet theory is a generalization over brane-gravity. 
So, it is reasonable to analyze this scenarion in in GB-model also. 
Modified Priedmann equations for Einstein- Gauss- Bonnet theory is 

given as 



4a 



(l-4(5At')cos/i(y) - 



(34) 



{gs is the string energy scale), energy density P[gb) and 



where a 

brane-tension X(^gb) is connected to energy parameter x as 

r2(l -4a//2^)3^i/2 



kl{p{GB) + \gb)) 



a 



sinhx 



where 



^5A^3B) = 2/.(3-4a/x^), 
(1 + Aaii^) 2 Stt 



2 _ -r -^^1^ , 2 _ 

Ate — hjA — 



1^ 



Ml 



and 



k^ 



Stt 
Ml' 



(35a) 

(356) 
(35c) 

(35rf) 



y^QB) is obtained using the condition P{gb) — H — [15, 16]in eq.(35a). 
Like RSII model, here also, matter field is given by the conformal 

scalar field (b used above. So, P{gb) and P(gb) are given by eqs.(3c,d). 
In GB-model, conservation equation looks like 



P{gb) + '^H{p{GB) +P{gb) + \gb)) = 0. 



(36) 
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Connecting cqs.(35a,b,c,d) and (36), it is obtained that 



H 



klkl^ii^ + m 
' 2[^kl + ^diklH^] 



(P(GS) +P{GB) + \gb)) 



(37) 



Eqs.(3c) and (33a) yield 
= ^2 



(38) 



At H{tjn) — 0, using eqs.(3b,c,d), (35b) and (38) in eq.(37), it is 
obtained that 

1 + |(u;^ - 1)0^] H{tJ = -f [2(1 - y ) (Ag^ - Xgb) + Xgb] ■ (39) 
Case 1. For quintessence conformal scalar i.e. a;^ = 1 



In this case, eq.(39) reduces to 

H{tm) = 



kjU 



(40) 



It shows that unless Xqb — 4Ag^, H{tm) does not vanish, so barring 
this particular situation, in GB-model, the scale factor given by eq.(5a) 
is valid for quintessence conformal scalar too. It is unlike the case of 
GR-modcl as well as RSll-model. 

Case 2. For phantom conformal scalar (p i.e. a;^ = —1 



For this case, eq.(39) is obtained as 

1 - ^<i>'it^)\Hitm) = -| [^(Agi, - Xgb) + Xgb 



(41) 



This equation shows that H{tm) does not exist when (f)'^{tm) — 3/A;|. 
Moreover H{tm) vanishes when Xqb — (8/5)Ag^ and (t)^{tm) 7^ 3/A;|. 
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Thus, the scale factor a{t) given by eq.(5a) is vahd for phantom confor- 
mal scalar like GR-model as well as RSII-model, barring the situation 
when (f>\tm) = S/kj and Xgb = (8/5)Ag^ 
Quantum bounce near t — tg 

If expansion stops at t = tm, H{tm) = 0, so eq.(35a) yields 

PGs(tm)=Agjj-AGB. (42) 

It corresponds to GR-model in case Xgb = A^b- Quantum correction 
to eq.(35a) modifies the equation (42) as 

PGB{tm) = Xq]^ — Xgb — PA{tm) 
= ^GB ~ ^GB + p{tm) 

(43) 

where p{tjn) is given by eq.(17) for GR-model. This equation shows 
that pGB(^m) > p{tm), when Agg > Agb, otherwise pGB(^m) < p{tm)- 

Further universe contracts following the scale factor given by eq.(21a) 
and bounces &X, t — ts- At t = ts, the energy density is obtained as 

PGB{ts)^X'^B-^GB + p{ts), (44) 

where p{ts) is given by eq.(24a). 
Conclusion 

In the cosmic picture explored above, it is found that quantum grav- 
ity makes drastic changes in the course of future accelerated universe 
driven by phantom dark energy in GR-model and RSII-model. But ,in 
case of GB-model, similar results are obtained for both quintessence 
dark energy and phantom dark energy. It is demonstrated above that, 
under quantum effects, universe expands upto a maximum when energy 
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density grows to p{tm), PRs{tm) , and pcBitm), given by eqs.(17), (32) 
and (43) in GR , RSII and GB- models respectively. It is found that, 
in RSII model,due to brane-tension, pRsitm) is lower than the same 
in GR model. In GB- model, pGB{tm) is lower than the same in GR 
model if brane-tension is higher than the critical brane-tension, but for 
brane-tension lower than the critical brane-tension pGB{tm) is higher 
than p{tm)- 

Later on, it contracts upto a minimum size Us at t — ts, when energy 
density rises to ps, given by eqs.(24a), (33) and (44) for GR , RSII and 
GB- models respectively. Like pRs{tm): PRs{ts) is also lower in RSII- 
model. In GB-model, pcsits) is lower or higher than p{ts) depending on 
the brane-tension like pGB{tm)- Finite energy density is obtained when 
a — (3 — 3 and 7 = 4 in the scale factor a{t) given by eq.(18a). It shows 
that universe will contract during the time period tm < t < tg and 
expand for t > ts exhibiting a bounce at t — ts- Thus, a possibility of 
re-birth to an exponentially expanding non-singular universe for t > tg 
is obtained. In the present cosmic picture, future singularity is not 
mild, but cosmic bounce at t — tg shows its complete avoidance. 
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